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Abstract 

We describe a computational approach to the verification of Maeda's conjecture 
for the Hecke operator T2 on the space of cusp forms of level one. We provide 
experimental evidence for all weights less than 12 000, as well as some applications 
of these results. The algorithm was implemented using the mathematical software 
Sage, and the code and resulting data were made freely available. 

P 1 Introduction 

Modular forms come in many different types. One of the most attractive aspects of 
the theory is that, despite the apparent variety of definitions and properties, there are 
some universal guiding principles (such as the Langlands program) that serve to unify 
and motivate this diversity. On the other hand, there are some special properties that 
seem to occur in isolation. One such instance is provided by a conjecture formulated by 
1— i Maeda, which indicates a behavior that seems to be specific to modular forms of level 

one on GL2. 

Before describing Maeda's conjecture in more detail, we review some basic definitions 
and properties of modular forms. For a thorough treatment of the background needed 
in this paper, the reader is invited to consult |Ste07 



Let k G 7L. A modular form of level 1 and weight /c is a holomorphic function 
C\j f:H — >C, where % = {z G C | Imz > 0}, 

£al satisfying 

• Modularity: for all z G H and all g = G SL 2 (Z), 
d ( az + b \ k 

• Holomorphicity at ico: a holomorphic function / satisfying the modularity condi- 
tion satisfies f(z + 1) = f(z) for all z G 7i, so it has a Fourier expansion 



00 

f(z) = a n q n , where we set q = e 2mz 
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We ask for / to be holomorphic at too, i.e. that a n = for all n < 0. 

We say that a modular form / is a cusp form if ao = 0. The cusp forms of weight 
k form a vector space Sjfc. These vector spaces are equipped with a family of Hecke 
operators T m (for m € N), whose effect on the Fourier expansion f(q) = ^a n q n of 
/ G Sfc is given by 



(T m /)(?) = E E d*" 1 *™/* k 

n=l W|gcd(m,n) 



The complex vector space 5^ has dimension 

- 1 if k = 2 
if fc ^ 2 




(mod 12), 
(mod 12). 



Let F denote the characteristic polynomial of the operator T2 acting on St , and let 
d = dimSfc. In the 1970s, Yoshitaka Maeda noticed that F is irreducible over Q for 
all k such that al < 12. In the 1990s, Lee-Hung |LH95 and Buzzard Buz96 studied 



these polynomials further and observed in a number of cases that the Galois group of F 
is the symmetric group &d- Shortly thereafter, Maeda made the following conjectural 
statement: 



Conjecture 1.1 (Maeda |HM97|). Letm > 1 and let F be the characteristic polynomial 
of the Hecke operator T m acting on . Then 

(1) the polynomial F is irreducible over Q; 

(2) the Galois group of the splitting field of F is the full symmetric group &d, where d 
is the dimension of . 

The conjecture has enjoyed constant attention over the last 15 years, with theoreti- 
cal as well as computational results. We summarize the computational verifications in 
Tabled) 



Source 



l g |LH 95| 
. |Buz96| 
HM97T 



Lee-Hung JLH95 
Buzzard Buz96 
Maeda 



CF99 



FJ02 



Conrey- Farmer 
Farmer-James 
Buzzard-Stein, Kleinerman 
Chu-Wee Lim 
present paper 



Lim05 



Kle04 



weights 



k < 62, k + 60 
k = 12£, £ prime, 2 < I < 19 
k < 468 

k < 500, k = (mod 4) 
k < 2 000 
k < 3 000 
k < 6 000 
k < 12 000 



Table 1: Summary of known cases of Maeda's conjecture for T2 

The theoretical results focus on whether the validity of the conjecture for a given 
operator T m can be used to deduce the conjecture for other operators T n . We state three 
such results, each giving a partial answer to this question. 
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Theorem 1.2 (Conrey- Farmer- Wallace CFW00| ). Let k be a positive even integer. 



Suppose there exists n > 2 such that the operator T n acting on Sk satisfies Maeda's 
conjecture. Then so does T p acting on Sk, for every prime p in the set of density 5/6 
defined by the conditions 

p ^ ±1 (mod 5) or p ^ ±1 (mod 7). 

Stated differently, this says that if Maeda's conjecture in weight k holds for one index 
n, then the density of primes for which the conjecture fails is at most 1/6. The next 
result considers only the irreducibility part of the conjecture, but it is stronger since it 
says that the density of primes for which the conjecture fails is zero. 

Theorem 1.3 (Baba-Murty |BM03|). Let k be a positive even integer. Suppose there 
exists a prime p such that the characteristic polynomial ofT p acting on Sk is irreducible 
over Q. Then there exists 5 > such that 

N 

< N prime I charpolyiTASk) is reducible} <C — x . 

(log N) 1+0 

Finally, Ahlgren gave a simple criterion for extending the validity of Maeda's conjec- 
ture from one index to another, and used it together with some computer work to prove 
the following result. 

Theorem 1.4 (Ahlgren |Ahl08|). Let k be such that d := dimS^ > 2. Suppose there 
exists n > 2 such that the operator T n acting on Sk satisfies Maeda 's conjecture. Then 

(1) T p acting on Sk satisfies Maeda's conjecture for all primes p < 4 000 000/ 

(2) T n acting on Sk satisfies Maeda's conjecture for all n < 10 000. 
We can now state our main result. 

Theorem 1.5. Let k < 12 000 and let 

n e{2, . . . , 10 000} U {p prime | 2 < p < 4 000 000} 

U {p prime \ p ^ ±1 (mod 5)} U {p prime \ p ^ ±1 (mod 7)}. 

Let F be the characteristic polynomial of the Hecke operator T n acting on the space Sk 
of cusp forms of weight k and level 1. Then F is irreducible over Q and the Galois group 
of its splitting field is the full symmetric group &d, where d is the dimension of the space 
Sk- 

Proof. The statement for T2 is the result of the computations described below. Given 
this, we deduce the result for the other T n by applying the results of Conrey-Farmer- 
Wallace and Ahlgren, as stated above. □ 

Our computational approach follows the "multimodular" method introduced by Buz- 



zard in |Buz96| and refined by Conrey- Farmer in CF99 . The main improvement is the 
use of random primes of moderate size, instead of going through primes consecutively 
until suitable ones are found. In Section [3] we describe the theoretical foundation of this 
approach, and we estimate the densities of the different types of primes we are looking 
for. This provides us with expected running times for our randomized algorithm, the 
Sage implementation of which we discuss in detail in Section |4j Finally, Section [5] gives 



some direct corollaries of Theorem 1.5 to some questions about modular forms of level 
one. 

We have made the code and data used to verify Theorem |1.5| available at 
http : / /bitbucket . org/ aghitza/maeda_data 



3 



2 Polynomial factorization and Frobenius elements 



Since this correspondence plays a central role in the theoretical underpinnings of our 
algorithm, we review some facts about factorization of polynomials over finite fields and 
Frobenius elements in Galois groups. These ideas go back all the way to the beginnings of 
algebraic number theory, appearing for instance in the work of Frobenius. A fascinating 
exposition of the mathematics and history of these ideas is given by Stevenhagen and 
Lenstra in [SL96 



We start with a bit of terminology. If r is a permutation on d letters, it can be 
decomposed into a product of disjoint cycles, uniquely up to permutation of the cycles. 
We say that r has cycle pattern d r [ ll d?, 1 ' 2 . . . dV 1 * if its decomposition contains exactly rrij 
cycles of length dj, for j = 1, . . . ,t. If H is a polynomial in F P [X], we say that H has 
factorization pattern d™ 1 d™ 2 . . . d™* if H has exactly mj irreducible factors of degree dj 
over F p . We recall that H is said to be separable if it has distinct roots over ¥ p . 



Lemma 2.1. Let F G 7L\X\ be monic, let p be a prime and let F p G F p LY] be the 
reduction of F modulo p. If F p is separable, then there exists an element r of the Galois 
group of F such that the cycle pattern of r is the same as the factorization pattern of 
F p . 

We sketch a proof of this classical result. 

Fix a prime p and consider the field automorphism a : ¥ p — > ¥ p given by o~(a) = a p . 
Since a fixes the subfield ¥ p , it permutes the roots of any polynomial H G ¥ P [X]. 
Moreover, Galois theory tells us that the cycle pattern of a (viewed as a permutation) 
is the same as the factorization pattern of H over ¥ p . 

We now take a monic polynomial F G %[X] and we let K/Q be its splitting field, 
Ok the ring of integers of K, and G the Galois group of K/Q. Let p be a prime in 
Ok over p. Suppose the reduction F p of F modulo p is a separable polynomial (in this 
case, we say that p is unramified in K/Q). Then there is a Frobenius element Frobp G G 
determined uniquely by the property 

Frobp(a) = a(a) (mod p) for all a G Ok- 

This implies that Frobp permutes the roots ai, . . . , ay G Ok of F in the exact same 
way as a permutes the roots in ¥ p of F p . We conclude that the cycle pattern of Frobp 
is the same as the factorization pattern of F p over F p . Therefore we can take r in the 



conclusion of Lemma 2.1 to be Frobp. 

Note that r is not uniquely determined by F and p, as the choice of a prime p of Ok 
above p matters. However, any two such r are conjugate in the Galois group. 



The following result follows easily from Lemma 2.1 and the fact that for any F G 
there are only finitely many primes p (namely the ones dividing the discriminant of F) 
for which F p is not separable. 

Theorem 2.2 (Frobenius). Let F G Z[X] be monic, let K/Q be the splitting field of F 
and let G be the Galois group of K/Q. Let degF = m\d\ + . . . + mtdt be a partition 
of deg F. The density of primes p for which F p has factorization pattern d™ 1 . . . d™ 1 is 
equal to 

#{<T G G | the cycle pattern of a is d™ 1 . . . d™'} 
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3 The basic lemma and density estimates 

Consider a monic polynomial F G %[X] of degree d. Given a prime p, we denote by 
F p £ Fp[X] the reduction modulo p of F. We say that the prime p is 

(1) of type I if F p is irreducible over ¥ p ; 

(2) of type II if F p factors over ¥ p into a product of distinct irreducible factors 

Fp = /o/i ■ • • fs 

with 

deg/ = 2 

deg/j odd for j = 1, . . . ,s; 

(3) of type III if F p factors over ¥ p into a product of distinct irreducible factors 

F p = /o/i ■ • • fs 

with deg/o > d/2 and prime. 



Remark 3.1. Hida and Maeda use a similar approach in Section 5 of HM97J, but replace 
primes of type III with primes of type IV, i.e. p such that F p = /o/i with /o, /i distinct 
and irreducible, and deg /o = 1. We will see below that primes of type III are significantly 
more common (and therefore better suited for our algorithm) than those of type IV. 

Remark 3.2. These types are not necessarily mutually exclusive: if d itself is prime, then 
a prime p of type I is clearly also of type III. 

Remark 3.3. In either of the three types, the conditions imply that the reduced polyno- 
mial Fp is separable: 

• If p is a prime of type I, then F p is irreducible, hence separable. 

• If p is a prime of type II or III, F p is a product of distinct irreducible factors. Each 
of the factors is then separable, and they cannot have any common roots, since 
otherwise they would have a nonconstant common factor and would therefore be 
reducible. Hence F p has distinct roots. 

Our computational approach to Maeda's conjecture is based on the following result, 
first proved in a special case in [Buz96| and then generalized in |CF99 . 



Lemma 3.4 (Buzzard, Conrey- Farmer) . Let F £ Z[X] be a monic polynomial of degree 
d. Suppose that F has primes of respective types I, II and III. Then F is irreducible over 
Q and its splitting field over Q has full Galois group &d- 

Proof. The fact that F is irreducible is immediate from the existence of a prime of type 
I. 

Let K/Q be the splitting field of F and let G be the Galois group of K/Q. Since F 
is irreducible, G is a transitive subgroup of &d- 



We also have a prime of type II. By Lemma 2.1 there exists t\ £ G whose decompo- 
sition into disjoint cycles contains exactly one even cycle (of length 2). Let a be the least 
common multiple of the lengths of the other cycles in t\, then rf £ G is a transposition. 
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Finally, there is a prime of type III. By Lemma 2.1 there exists T2 G G whose 
decomposition into disjoint cycles contains one cycle of prime length p > d/2. Therefore 
the other cycles have lengths that are coprime to p; letting b denote the least common 
multiple of these lengths, we find that r\ G G is a p-cycle. 

We now use the existence of these elements of G to conclude that G = &d- For 
i,j G S = {1, . . . ,d}, write i ~ j if i = j or if the transposition (i j) is in G. This is 
an equivalence relation on S. Since G is transitive, each equivalence class has the same 
number n of elements and it follows that n \ d = #S. Note that n > 1 since G contains 
at least one transposition, namely rf . Let T be the subset of S permuted by r|, and let 
Gt be the subgroup of G fixing S \ T. Define an equivalence relation on T by % ~ j if 
i = j or if the transposition (i j) G Gy. As before, each equivalence class has the same 
number m of elements and m \ p = #T. Since n > 1, we have m > 1, so m = p since p is 
prime. But n > m because Gt C G. Thus n > d/2, so n = d. This implies G = &d- D 

Our algorithm will consist of picking random primes and checking whether they are 
of type I, II or III for the characteristic polynomial of the Hecke operator T%. According 



to Theorem 2.2 , it is therefore important to estimate the number of permutations having 
certain types of cycle patterns. For a fixed pattern, the following well-known result (see, 
for instance, Proposition 1.3.2 of Sta97| ) gives an exact expression for the number of 



permutations. 

Lemma 3.5. Let an element a of &d have cycle type d 7 ^ L1 d™ 2 ...d™ 1 , where rrii is the 
number of times a cycle of length d{ appears in the cycle decomposition of a. (Note: 
midi+m2d2+- ■ .+mtdt = d.) The number of elements of&d of cycle type d^d™ 2 ■ ■ ■ d™ 1 
is equal to 

Cid™ 1 ^ 2 ... dT f } 



n U C™;' 



Proposition 3.6. The density of primes of type I is 

D M = \. 

Proof. Primes of type I correspond to d-cycles in &d- Each such cycle can be written 
uniquely as a sequence 1, a±, . . . , a^-i, where a%, . . . , ad-i G {2, . . . , d} can appear in any 
order. Therefore there are (d— 1)! d-cycles, and by Theorem |2.2| the density of primes 
of type I is 

(<*-!)! _ 1 
d\ d' 

□ 

In order to state our result on primes of type II, recall that for n G Z>o odd, the 
double factorial n\\ of n is the product of all the odd positive integers less than or equal 
to n. 

Proposition 3.7. Let d > 2 and let d be the largest even integer such that d < d. The 
density of primes of type II is given by 

V ; 2(d-2)\ 

and satisfies the inequality 

D " {i) > ITS' 
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Proof. Primes of type II correspond to elements in <3^ containing a 2-cycle and no other 
even cycles. There are Q 2-cycles in &d] fixing a 2-cycle, we need the number ff(d — 2) 
of elements of odd order in &d-2- We have 



d\\2J K ' 2(d-2)!' 



The sequence (£?(n) | n E N) appears in nature in several guises, see |OE IS|. The 



recurrence formulas that appear there and in Chapter IV of |Rio58| easily give the exact 
expression 

(n — 1)!!, if n is even 
(n — 2)\\n, if n is odd, 



which immediately provides us with the exact expression for Djj(d) in the statement. 

It remains to establish the lower bound. Write d = 2c for some c E Z (recall that d 
is the largest even integer less than or equal to d). Then 

m 2) _ [(2c-3)!!F _ (2c-3)l 

[ ' ~ 2(2c-2)! ~ 22«=-3(2c-2)[(c-2)!]2' {6A > 

We use the following bounds on the factorial, which can be thought of as an effective 



version of Stirling's approximation and were obtained by Robbins (see |Rob55 and 
Section II.9 in [Fel68] ): 



Then by using the lower bound for the numerator and the upper bound for the denom- 



inator on the right hand side of Equation (3.1), we obtain 



, , 2c - 3 1 i i /9\ 2 \/2 1 

D H (d) > e 2Mc-2) + l 6(c-2) > — 

' 2c-2 2^A(c-2) \WJ V^Vd 
where we use the elementary inequalities (valid for c > 6): 



g 24(c-2) + l 6(c-2) > p~2A c > 

2c - 3 9 
2c- 2 ~ 10" 

Since 

9 \ 2 V2 1 

10 J > 4' 



this gives us the desired lower bound for d > 12, and the remaining cases 2 < d < 12 are 
easily checked. □ 

Proposition 3.8. The density of primes of type III is 

d/2<£<d, I prime 

If d > 2, then 

D in (d) > i 
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Proof. Fix a prime £ such that d/2 < I < d. According to Theorem 2.2 we need to 
count the number of elements of 6^ that contain an ^-cycle. Choosing the ^-cycle itself 
involves the (fy ways of picking its constituents, which can then be rearranged within 
the cycle in {£ — 1)! ways. It remains to take into account the number of permutations 
of the remaining d — £ symbols, so overall we have 

A (i - iy.(d - i)\ = | 

elements of 6^ containing an ^-cycle, which gives the stated density. 

The inequality given in the statement follows from Bertrand's postulate (proved by 
Chebyshev), which says that for any integer n > 1 there is at least one prime £ such that 
n < £ < 2n. □ 

We can get a much better lower bound on the density Dm by using some recent 
results of Dusart on explicit estimates for sums over primes. 

Theorem 3.9 (Dusart, Theorem 6.10 in [DuslO] ). Let B & 0.26149 denote the Meissel- 
Mertens constant. For all x > 1 we have 

loglogx + £-( , 4 , ) < (3.2) 

Vl01og 2 x 151og 3 xy ^ x p 

We will also need an upper bound on the sum of the reciprocals of primes up to x, 
but Dusart's upper bound only holds for x > 10 372. For our purposes, the following 
weaker result is sufficient: for all x > 1 we have 

V- <loglogx + £ + -^. (3.3) 
^ p log z X 



(This inequality can be found in Theorem 8.8.5 of BS96| .) 
Proposition 3.10. If d > 10, then 

Dni(d) > ' 



3 log d 



Proof. We put together inequalities (3.2) and (3.3) to get 

d 1 4 

Dm (d) > log log d - log log ; 



2 10 log 2 d 15 log 3 c2 log 2 f 



2 

We write 



log log d — log log - = log ( 1 + 



2 ° V log d - log 2 
and use the inequality 



2 3 4 

X X X 



logfl + x) > x 1 for all 1 < x < 1 

&v 1 ~ 2 3 4 
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to get that for all d > 4 

D m(d) >- -1— — 

log a — log 2 



log 2 



log 2 2 11 



1 



1 



> 



log d 



0.693 - 1.341 
- 0.058 



2 10 y log d - log 2 
4 _ log 3 2\ 1 _ log 4 2 

15 3/ (log d- log 2) 2 ~~ 4 (log d- log 2) 3 
1 1 



log d — log 2 
1 



0.156 



(log d- log 2) 2 



(log<i — log2) s 



If d > 94, then the expression in the brackets is bigger than 1/3, and we get the desired 
inequality. We check that it holds for the remaining cases 10 < d < 94 by computation. 

□ 



For completeness, we treat the case of primes of type IV, as defined in Remark 3.1 
Proposition 3.11. Let d > 1. The density of primes of type IV is 

Proof. We need to count the number of (d — l)-cycles in &d- There are d choices for the 
letter that is fixed, and (d— 2)! choices for permuting the other letters appropriately, 
therefore the density of primes of type IV is 



d{d-2)\ 
d\ 



1 

d- 1" 



□ 



4 Implementation and results 



Our approach is a randomized version of the algorithm from |CF99 , based on the results 



introduced in the previous section. We implemented this algorithm using the mathe- 



matical software Sage, see |Sagc 



Here is a description of the main steps used to verify Maeda's conjecture for a fixed 
weight k; in those cases where a major step is delegated to a component of Sage (rather 
than using native Sage code), we mention the relevant component. 

(1) Compute the Victor Miller basis B for Sk up to precision 2(d + 2), where d is the 
dimension of The Sage implementation of this basis uses |FLINT] polynomials 
as the internal data structure. 

(2) Compute the matrix M of the Hecke operator T2 with respect to the basis B - this 
is very efficient since the basis B is echelonized. 

(3) Pick a random prime p < 2 20 , uniformly over this range. (This choice of upper 
bound gives a large enough range so that it is likely to contain primes of type we 
are looking for, but not so large that the arithmetic over ¥ p gets too expensive.) 

(4) Reduce M modulo p and compute the characteristic polynomial F p E F P [X]. The 



characteristic polynomial is computed by the Linbox| library. 
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(5) Is F p irreducible? If so, p is a prime of type I. The irreducibility test uses FLINT . 



(6) Factor F p over ¥ p and use this factorization to decide whether p is a prime of type 
II or III. The factorization is done by [FLINT] . 

(7) Repeat from step ^ until we have found at least one prime of each type. 



According to Propositions 3.6 3.7 and 3.10 we expect to look on average at d primes 
before we find one of type I, at A^fd primes to find one of type II, and at 3 log d primes 
to find one of type III. 

The actual performance of this algorithm (as well as a comparison to the consecutive 
version of the algorithm, used in [ CF99] ) is illustrated in Figure [T] Some care needs to 



be taken in interpreting the graphs: 

• There is no difference in running times for Steps ([I]) and ([2]), which are common 
between the two algorithms. 

• As the weight increases, the major component of the running time is finding a prime 
of type I. Therefore, even though the randomized algorithm does much better at 
finding primes of types II and III, this advantage has only a minor impact on the 
overall running time. 

• In the range illustrated in the graphs (i.e. weights less than 2 000), the randomized 
algorithm required on average one third of the number of primes needed by the 
consecutive algorithm. However, some of this is counteracted by the fact that the 
consecutive algorithm works with much smaller primes, which are faster to test. 

• Overall, for weights less than 2 000, the randomized algorithm was about twice as 
fast as the consecutive one. 

For weight k = 12 000, the entire verification took about 46 hours. The majority 
(93%) of the time was spent looking for a prime of type I; this required testing 1 783 
primes, and each test took about 87 seconds. The computation of the Victor Miller basis 
took about 2.2 hours, and the computation of the characteristic polynomial of T<i took 
about 1 hour. 



5 Some applications 



We record some immediate consequences of Theorem 1.5 



5.1 Non- vanishing of L-functions 

A modular form is called an eigenform if it is an eigenvector for all the Hecke operators 
T n . The L-function associated to an eigenform / = Y2ri=i a nQ n of weight k is given by 

oo 

«/,.)- 

n=l 

If k = 2 (mod 4), the functional equation of L implies that L(f, k/2) = 0. It is believed 
that if k = (mod 4), then L(f,k/2) ^ 0. The following result follows immediately 
from work of Conrey- Farmer: 

Corollary 5.1 (see Theorem 1 in |CF99] ). Suppose k = (mod 4) and k < 12 000. 

Then L(f, k/2) ^ for any cuspidal eigenform f of level 1 and weight k. 
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5.2 Base change for totally real fields 



It is in the context of this work of Hida and Maeda that Maeda's conjecture was formu- 
lated. We content ourselves with giving a general description of this application, and we 



refer the interested reader to HM97 for details. 

Let / E Sk be a Hecke eigenform. For each prime p, there is a p-adic Galois repre- 
sentation 

p: Gal (Q/Q) ^GL 2 (Q p ). 

There is an Artin L-function L(p, s) attached to p, and the relation between p and / 
can be summarized by 

L(p,s) = L(f,s). 

Now let E be a number field. There is a purely algebraic notion of a cohomological 
eigenform / on GL^Ag), where A^ is the ring of adeles of E. We say that / is a base 
change of f to E if 

L{f,s) = L(p E ,s), 
where pe ■ G&\(Q/E) — ► GL2(Q„) is the restriction of p to E. 



The work of Hida and Maeda, together with Theorem 1.5 implies that for k < 12 000 
and a totally real field E satisfying some ramification conditions, any eigenform f £ Sk 
has a base change to E. 

5.3 Eigenforms divisible by eigenforms 

It is easy to see from the definition of a modular form that if f\ and fi are modular 
forms of respective weights fci and &2, then the product /1/2 is a modular form of weight 
k\ + &2- In other words, modular forms of all weights put together form a graded algebra 



M = 0M fc . 



A natural question is whether the product of eigenforms can be an eigenform. This 
will clearly happen for small weights (for instance, when the product lives in a one- 
dimensional space of cusp forms). Since the Hecke operators do not act on the entire 
algebra M of modular forms (they act differently on the graded pieces Mfc), it seems 
reasonable that the one- dimensional coincidences are the only situation in which a prod- 
uct of eigenforms is an eigenform. Such questions have been studied by several authors, 



with the latest results appearing in a recent paper by Beyerl-James-Xue |BJX11 . They 



consider the more general question of divisibility of an eigenform by another eigenform, 
i.e. relations of the form h = fg where /, g, h are modular forms and /, h are eigen- 
forms. The relation with Maeda's conjecture is discussed in Section 6 of |BJX11 , and 
Theorem 1 1 . 5| implies the following result. 



Corollary 5.2. Let h be a cuspidal eigenform of weight < 12 000, and let f be an 

eigenform (which could be cuspidal or Eisenstein). Then h = fg for some modular form 
g E Mfc with k > 2 if and only if we are in one of the cases listed in Table [S| 



5.4 Distinguishing Hecke eigenforms 

How many initial Fourier coefficients are necessary to completely determine a Hecke 
eigenform? Theorem 1 in |Ghill says that a 2 , 03 and 04 are sufficient, but our compu- 
tational verification of Maeda's conjecture gives a stronger result: 
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weight of / 


weight of g modulo 12 


nature of / 


4 


0,4,6, 10 


Eisenstein 


6 


0,4,8 


Eisenstein 


8 


0,6 


Eisenstein 


10 


0,4 


Eisenstein 


12 


0,2,4,6,8, 10 


cuspidal 


14 





Eisenstein 


16 


0,4,6,10 


cuspidal 


18 


0,4,8 


cuspidal 


20 


0,6 


cuspidal 


22 


0,4 


cuspidal 


26 





cuspidal 



Table 2: The only cases in which a cuspidal eigenform of weight < 12 000 can be factored 



into h = fg with / an eigenform, see Corollary 5.2 



Corollary 5.3 (see Theorem 6 in |Ghill| ). Comparing the Fourier coefficient 02 is 
sufficient to distinguish any two cuspidal eigenforms of level 1 and weights < 10 000. 
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Figure 1 : Histograms illustrating the number of primes tested before finding a prime of 
type I, II, respectively III, in weights up to 2 000. In each graph, the numbers on the 
x-axis represent the ratio N/E of the actual number of primes tested over the expected 
number of primes (coming from the densities described in Section [3]). The y- value 
represents the number of weights featuring (a small neighborhood of) that particular 
ratio N/E. The blue continuous line corresponds to our randomized algorithm, while the 
red dotted line corresponds to the consecutive algorithm from |CF99| . As an example: in 
the top graph, the global maximum on the continuous line is at (0.1, 101), meaning that 
for 101 weights, the number of candidates for a prime of type I tested in the randomized 
algorithm was about 1/10 of the expected number of primes. 
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